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We study the origins of period doubling and chaos in the 
driven series resistor-inductor-varactor diode (RLD) nonlin- 
ear resonant circuit. We find that resonators driven at fre- 
quencies much higher than the diode reverse recovery rate 
do not show period doubling. Models of chaos based on the 
nonlinear capacitance of the varactor diode display a reverse- 
recovery- like effect, and this effect strongly resembles reverse 
recovery of real diodes. We find for the first time that in ad- 
dition to the known dependence of the reverse recovery time 
on past current maxima, there are also important nonlinear 
dependencies on pulse frequency, duty-cycle, and DC voltage 
bias. Similar nonlinearities are present in the nonlinear ca- 
pacitance models of these diodes. We conclude that a history- 
dependent and nonlinear reverse recovery time is an essential 
ingredient for chaotic behavior of this circuit, and demon- 
strate for the first time that all major competing models have 
this effect, either explicitly or implicitly. Besides unifying the 
two major models of RLD chaos, our work reveals that the 
nonlinearities of the reverse recovery time must be included 
for a complete understanding of period doubling and chaos in 
this circuit. 



I. INTRODUCTION 

The nonlinear dynamics of the p-n junction has been 
a subject of intense interest since the dawn of chaos. 1-8 
Experimental work has focused on the damped driven 
nonlinear oscillator formed by a resistor, inductor, and 
varactor diode (RLD circuit) connected in series (see 
Fig. 1(a)). This is the simplest passive circuit that dis- 
plays period doubling and chaos as a function of sinu- 
soidal driving amplitude and frequency. The dynamics 
of this circuit can, in principle, be understood through 
the details of the hole and electron charge distributions in 
the vicinity of the p-n junction. 9 However, many simpli- 
fied models of p-n junction charge dynamics have been 
put forward in terms of nonlinear lumped-element ap- 
proximations for the varactor (Fig. 1(b)). One reason 
for this approach was the interest in constructing low- 
dimensional maps to describe the "universal" nonlinear 
dynamics experimentally observed in the circuit. 1,2,4,6 
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FIG. 1. Schematic diagrams for the (a) driven RLD cir- 
cuit, and (b) diode models. The nonlinear capacitance and 
nonlinear resistance (NLCR) implementation is shown in the 
upper part of (b), while a representation of the Hunt, etal. 
model is shown in the lower part. In the latter case, the diode 
is replaced with a battery while it is conducting, and a fixed 
linear capacitor at all other times. 
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There are three general types of nonlincarity present 
in the driven RLD circuit. The first is the nonlinear 
current-voltage characteristic of the diode, thought by 
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most researchers to be unimportant for chaos, because 
its absence from models results in only small changes in 
the calculated bifurcation diagrams. 5 ' 10 A second con- 
tribution comes from the large and exponentially nonlin- 
ear forward-bias capacitance associated with the junction 
diffusion capacitance. 5 In this case, the diode is modeled 
as a parallel combination of a nonlinear resistor and non- 
linear capacitor (Fig. 1 (b)). It has been proposed that 
period doubling will occur when the capacitance reaches 
approximately four times its zero bias value because the 
resonant frequency of the RLD circuit will drop to half its 
low-amplitude value. 7 This is thought to be the first step 
in a period doubling cascade to chaos. A third contribu- 
tion to nonlinearity comes from the finite time-scale diffu- 
sive dynamics of charge in the p-n junction and the asso- 
ciated "memory" of previous forward-current maxima. 4 
After the diode has been forward-biased and switched 
off, it takes some time for the minority carriers to diffuse 
back across the junction, allowing the diode to conduct 
for a period known as the reverse-recovery time, trr. 
The "memory" of previous forward bias currents built 
into this non-equilibrium charge distribution has been 
proposed as the main source of chaos in the driven RLD 
circuit. 4,8 

Lumped element models of the driven RLD circuit fall 
into two distinct classes: those relying on the nonlin- 
ear capacitance and ignoring explicit history effects (Fig. 
1(b) NLCR), and those explicitly based on history effects 
in the junction and ignoring the nonlinear capacitance 
(Fig. 1(b) Hunt). Both kinds of model attempt to iden- 
tify the essential physics that is responsible for period 
doubling and chaos. It is remarkable that both classes 
of models reproduce intricate detail of the return maps 
and bifurcation diagrams of the driven RLD circuit. His- 
torically, the two classes of models were assumed to be 
mutually exclusive and contradictory in their approach to 
describing chaotic dynamics in the driven RLD circuit. 11 
One purpose of this paper is to partially reconcile these 
two classes of models, and to show that they both include 
some essential physics of this rich nonlinear problem. We 
also demonstrate that both models lack important effects 
that have been ignored up to this point. These effects 
are measured here for the first time in the context of 
nonlinear dynamics, and are essential for a complete un- 
derstanding of the behavior of this circuit. 

In this paper we present new measurements and simu- 
lations of reverse-recovery times in a variety of diodes 
and resonant circuits. Section II discusses the many 
time scales of the driven RLD problem and illustrates 
the importance of the reverse recovery time in the non- 
linear dynamics of the circuit. Section III demonstrates 
that the nonlinear capacitance models are endowed with 
a history-dependent reverse-recovery-like effect. Section 
IV is a summary of our experimental results of the non- 
linear dependence of the reverse recovery time on circuit 
parameters, and ends with a dramatic illustration of how 
existing models fail to explain key data. Sections V and 
VI review our contributions to a deeper understanding of 



the RLD circuit and demonstrate that a simple nonlin- 
ear capacitor model shows all of the essential features of 
reverse-recovery-like behavior. 

Experiments were carried out in a manner similar to 
that described in Ref. 12 . The circuits were shielded in- 
side a metallic box and excited through 50 f2 transmis- 
sion lines by either an HP 33 120 A synthesized source 
(up to 15 MHz) or an HP 83620B microwave synthesizer 
(above 10 MHz). The response of the circuits was mea- 
sured with either a Tektronix TDS3052 oscilloscope, or 
a Tektronix 494P spectrum analyzer. Low capacitance 
and high impedance probes were used to measure volt- 
ages in the circuits. SPICE simulations showed that 
the probes had a minimal impact on the circuit wave- 
forms, whereas it was found that low impedance capac- 
itive probes strongly perturbed the circuits. Bifurca- 
tion diagrams were constructed by collecting single-shot 
traces from the oscilloscope and recording the data with 
a Lab View program on a PC. 

II. TIME SCALES 

The history of research on the driven RLD circuit is 
long and varied. Most papers on this subject study only 
one diode and choose a circuit with a resonant frequency 
in a convenient range (typically at audio frequencies). 
The experiments have been carried out in widely dis- 
parate regions of the parameter space available for the 
circuit (the parameters include driving frequency, ampli- 
tude, and duty-cycle, R and L values, minority carrier 
lifetime, dc voltage bias, etc.). 13-17 Only recently has 
effort been made to assemble a unifying picture of the 
origins and generality of period doubling and chaos in 
this important circuit. 10 Our purpose in this section is 
to measure and model different diode circuits in a wide 
range of parameter space to identify the important and 
essential physics responsible for the nonlinear dynamics 
of the driven RLD circuit as a step in creating such a 
unified picture. 

It is important to first understand the four major 
time/frequency scales involved in this problem, because 
the dynamics of the circuit depends critically on the rel- 
ative values of these time scales. The first is the low- 
amplitude resonant frequency of the passive RLD circuit: 
fo = l/p^LO,) 1 / 2 ], where Cj is the zero-bias junction 
capacitance of the varactor diode. It is well established 
that complex dynamics are most easily observed for driv- 
ing frequencies f in the vicinity of this resonant frequency. 
Another important time scale in the diode is the reverse 
recovery time: trr. This is a measure of the time re- 
quired for minority carriers injected across the junction 
in the diode to move back to the other side of the junction 
(or recombine) once the driving force has been removed. 
One effect of this finite time scale is "reverse recovery" 
in which a forward-biased diode will not shut off imme- 
diately when the driving signal is reversed, but will con- 
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tinue to conduct for a time on the order of trr. Su, 
Rollins and Hunt proposed that the diode actually acts 
like a battery during this time. In addition, they pro- 
posed that Tun retained a memory of the two or three 
previous maximum forward current values through the 
diode. 8 

There are several RC time constants in the problem, 
associated with the junction capacitance and the internal 
resistance of the diode, and between the junction capac- 
itance and the external resistor, R (note that the para- 
sitic resistance of the inductor is much less than R for the 
circuits considered here). These time scales are usually 
quite short (~ ns), but can be important when the capac- 
itance is large. Finally, the L/R time scale is generally 
long and does not directly influence the diode dynam- 
ics. The interplay of these time scales, along with the 
period of the sinusoidal driving signal, yield a very rich 
parameteric landscape for nonlinear dynamics. 

It should be mentioned that three of these time scales 
are strongly nonlinear. The resonant frequency fo and 
the RC time constants are a nonlinear function of driving 
amplitude and DC voltage bias on the varactor junction, 
and we show for the first time below, the reverse recovery 
time is a nonlinear function of drive frequency, ampli- 
tude, duty-cycle, DC bias, and load resistance. 9 Most of 
these nonlinearities have not been investigated in the past 
with respect to the resulting nonlinear dynamics. These 
nonlinearities greatly enrich the dynamics and challenge 
our understanding of the essential physics responsible for 
nonlinear dynamics in the simple driven RLD circuit. 

To illustrate the interplay of time scales in this prob- 
lem, consider the sinusoidally driven resonant RLD cir- 
cuit shown in Fig. 1(a). The low-amplitude resonant fre- 
quency of the circuit is determined by the diode junction 
capacitance Cj as well as the lumped element inductor, L. 
The reverse recovery time trr can be measured by meth- 
ods discussed below. We studied four different diodes 
with very different values of Cj and trr. In all cases we 
find the following. First, when the low-amplitude reso- 
nant frequency of the circuit fo is on the order of 1/trr, 
period doubling and chaos are found for a substantial 
driving frequency range f between about 0.2 fo and 2 
fo (see Table I). Similar results are obtained when the 
inductance is changed to make the resonant frequency 
equal to about one tenth of 1/trr (in agreement with 
Ref. 10 ). However, when the inductance is changed to 
make fo about ten times larger than 1/trr, there is a 
substantial change in the circuit behavior (Table I). One 
of the diode circuits shows only period doubling and no 
chaos, while another displays nothing but period-one be- 
havior. The other diode circuits show period doubling 
and chaos, but over a more restricted range of ampli- 
tude and driving frequency parameter space. Finally, 
when the circuits are modified to have fo ~ 100/r^, we 
do not observe period doubling or chaos at all over the 
equivalent range of parameters. 

When the diode is driven at frequencies much higher 
than the reverse recovery time, the charge carriers near 



the p-n junction do not have a chance to influence the cir- 
cuit behavior or contribute additional nonlinear behavior. 
Clearly the reverse recovery time plays an important role 
in the nonlinear dynamics of this circuit. However, very 
little prior experimental work on RLD dynamics has ex- 
plored the interplay of all the time scales in this problem. 



III. THE STORAGE (r s ) AND REVERSE 
RECOVERY (t rr ) TIME SCALES 

Hunt and collaborators claimed that the reverse re- 
covery time is the key feature that generates chaos in the 
driven RLD circuit because it is a function of the peak 
current flow in previous cycles, thus introducing a mem- 
ory effect into the dynamics. In the paper by Rollins and 
Hunt 4 a simple model of reverse recovery was introduced 
for RLD dynamics, and turned into a ID noninvertible 
map. Although this gave qualitative agreement with ex- 
periment, it was criticized for being inconsistent with the 
measured attractor. 5 Later, Su, Rollins and Hunt 8 in- 
cluded three or more previous cycles of the forward bias 
current to determine their model trr. They also noted 
that the recovery time depends on the junction charac- 
teristics and the external series resistance. We show here 
that this model is a good starting point, but fails to in- 
clude essential nonlinearities. 
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FIG. 2. (a) Current flowing through the diode during and 
after the pulse excitation shown in (b), demonstrating the 
reverse recovery effect for a 1N4007 diode (solid line) , and the 
discharge effect on a piecewise linear capacitor (PLC) model 
in place of the diode (dashed line). For the diode circuit 
R=25 O, while for the PLC 6 R = 60 O, Ci = 0.1 fiF, C 2 = 
400 pF, and E = 0.1 V (see Fig. 3(b)). These values for the 
PLC model are used throughout the paper, (b) Shows the 
drive voltage applied to both circuits and the inset shows the 
circuit schematic. 

Figure 2 illustrates the reverse recovery effect in a 
diode. The circuit shown in the inset of Fig. 2(b) was 
employed to drive the diode with a single square pulse 
input. After the driving signal was removed, the diode 
continued to conduct for a time scale defined as the re- 
verse recovery time, trr. This time consists of two con- 
tributions: the storage time r s , and the transition time 
r t , with trr = r s + r t . 9 We shall present results below 
(all data in Figs. 3-7 are taken with the circuit shown 
in Fig. 2(b)) mainly on the storage time because, in gen- 
eral, this time scale dominates the reverse recovery time 
and its nonlinear behavior. 

Figure 3 illustrates the history dependence of the stor- 
age and reverse recovery times. In this case we compare 
single pulse recovery data like that shown in Fig. 2 to 
the same diode subjected to two consecutive pulses (Fig. 
3(a)). The storage time is longer for the diode that has 
been subjected to two pulses, mainly because it was not 
finished recovering from the first pulse before the second 
arrives. This simple history dependence of the recovery 
time is the basis for the Hunt and Rollins treatment of 
complex RLD dynamics. 4,18 
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FIG. 3. Illustration of reverse-recovery history dependence 
through measurements of the voltage drop on the device 
(diode or PLC) after either one or two pulse excitation (taller 
curves) using the circuit in the inset of Fig. 2(b). The lighter 
curves show double pulse exposure and the darker curves show 
a single pulse exposure, (a) Diode 1N4007 and R=l kf2, (b) 
PLC model (implemented in SPICE) using the same circuit 
parameters as in Fig. 2. The PLC voltage-charge curve is 
shown in the inset. Note that in both cases the device requires 
a longer time to recover after the second pulse compared to 
just one pulse. 

All other models of RLD dynamics do not consider an 
explicit history dependence of the junction dynamics of 
the diode. They are based on the nonlinear behavior of 
the capacitor used to approximate the reactive electrical 
behavior of the diode (Fig. 1(b), NLCR). As we shall 
show here, in fact all of these other models do indeed have 
a history-dependence of a reverse-recovery-like time that 
naturally comes from the charge storage on the capacitor 
part of the varactor diode model. 

A. Piece-wise Linear Capacitor (PLC) Model 

To illustrate our point, we shall choose the simplest 
model of the nonlinear capacitance in the varactor. In 
this case, one can use a piece-wise linear approximation 
for the capacitance of the varactor, with just two lin- 
ear pieces. This model for the varactor capacitance was 
considered by several groups 6 ' 19 ' 20 and shown to produce 
period doubling and chaos when part of an R-L-PLC se- 
ries circuit driven by a sinusoidal signal. The inset to 
Fig. 3(b) shows that the PLC switches at potential E 
from a large value Ci for higher bias, to a small value 
C2 for lower bias. We use this model to illustrate that it 
has implicitly included a history-dependent storage time. 
Fig. 2(a) shows a SPICE simulation of the diode current 
(dashed line) using the PLC model for the varactor in the 
circuit shown in the inset of Fig. 2(b), where the diode 
is replaced with the PLC. It clearly displays a "storage 
time" and a "transition time", just like a real diode. 21 
The "storage time" (r s ) in this case corresponds to the 
RC time constant of the capacitor discharging into the 
external resistance of the circuit while it has a large ca- 
pacitance Ci. The "transition" part of the signal occurs 
when the device switches to its smaller capacitance value 
C2, producing a much shorter RC time constant that dis- 
charges the capacitor very quickly. Thus the discharge 
behavior of the capacitor mimics the effects of minority 
charge carrier diffusion across the junction. As we shall 
see, the details of this discharge behavior are remarkably 
similar to the diffusive dynamics of the minority charges 
in real diodes. 

We have found that a variety of models based on non- 
linear capacitance of the diode show a history dependence 
of the RC storage and "reverse recovery" times. Fig. 3(b) 
illustrates the history-dependence of the RC discharge 
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time in the PLC model. All models of RLD dynam- 
ics that we have investigated have a history-dependent 
reverse recovery effect coming from the nonlinear capac- 
itance. Our conclusion is that a finite reverse recovery 
time is indeed an important concept for understanding 
chaos in this circuit. Besides making this point in this 
paper, we also wish to explore the nonlinear and history- 
dependent behavior of this recovery time and how it in- 
fluences the generic behavior of the driven RLD circuit. 



IV. NONLINEAR BEHAVIOR OF THE STORAGE 

TIME 

The storage time and its history dependence are clearly 
important ingredients in any model of nonlinear dynam- 
ics of the driven RLD circuit. However, it is not widely 
appreciated that the storage time is itself a strongly 
nonlinear function of many parameters in the problem. 
These parameters include the amplitude of the forward 
bias current, the circuit resistance R, the driving fre- 
quency, the duty cycle of the drive, and the DC volt- 
age bias on the junction. The purpose of this section 
is to present experimental evidence and models of these 
nonlinearities and to discuss their significance for the oc- 
currence of period doubling and chaos in the parameter 
space of these circuits. 

First is the forward current amplitude dependence of 
the storage time r s . Hunt and Rollins pointed out that 
the reverse recovery time of a diode increases as the am- 
plitude of the forward bias current through the diode 
increases. The approximate expression used to describe 
this effect was; 4 ' 18 



trr = T m (l - exp(-\I max \/I c )), 



(1) 



where l max is the most recent maximum forward current, 
T m and I c are lifetime and current scales particular to 
each diode. A more exact expression for the storage time 
t s can be derived from analysis of charge dynamics in the 
p-n junction; 22 
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where Ir is the reverse current through the diode during 
the storage phase. All of these forms basically say that 
the storage and reverse recovery times will increase as 
the amplitude of the forward current increases. This is 
because more minority charge will be pushed across the 
junction by the larger current and, because of the diffu- 
sive dynamics of the carriers, there will be a longer delay 
in getting back to the equilibrium charge configuration. 



Data and PSpice simulations for the forward current 
dependence of the storage time in two different diodes 
is shown in Fig. 4. (Note that PSpice simulations use 
standardized semi-empirical models of the diode behav- 
ior, while our SPICE simulations of the PLC use the 
model shown in the inset of Fig. 3(b) with parameters 
given in the caption of Fig. 2.) The circuit in the in- 
set of Fig. 2(b) was used to perform these measurements 
and simulations. For the data in Figs. 4 and 5 the driv- 
ing signal is a periodic square wave with an equal time 
spent at the positive and negative amplitude voltage val- 
ues (50% duty cycle). The storage time was measured as 
a function of the driving signal amplitude, proportional 
to the forward current through the diode. The data and 
simulations all show an increase of the storage time with 
increased forward bias current. In some cases the storage 
time saturates, consistent with the approximate expres- 
sion Eq. (1), and in other cases it continues to rise slowly, 
consistent with the more sophisticated treatment in Eqs. 
(2) and (3). 
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FIG. 4. Squared wave amplitude dependences: a) Ampli- 
tude dependence of the storage time for NTE610 diode with 
(A) R=l kQ experiment, (B) R=l kQ. PSpice model, (C) 
R=25 Q PSpice, and (D) R=25 fi experiment, b) Ampli- 
tude dependence of the storage time for 1N4007 diode with 
(E) R=l kfl PSpice, (F) R=25 fi PSpice, (G) R=l kfi ex- 
periment, and (H) R=25 £7 experiment. Measurements and 
simulations were done using the inset circuit shown in Fig. 
2(b) where the driving voltage was a square wave shape (with 
no DC bias) of frequency 20 kHz (as the recovery time is 
constant for low frequencies). The peak-to-peak voltage was 
twice the amplitude. 

It has also been noted that the value of the resistance 
R in the circuit of Fig. 2(b) (and therefore Fig. 1(a)) will 
also influence the storage and reverse recovery times. 8 ' 18 
This is illustrated in Fig. 4. When the load resistance 
R is changed from 25 to 1 kfl, the storage time in- 
creases for both diodes (NTE610 and 1N4007) in both 
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experiment and simulation. This is expected because the 
reverse recovery process involves charge transport across 
the junction and through the rest of the circuit. A larger 
external resistance in the circuit will imply a smaller re- 
verse recovery current, thus extending the time for the 
minority charge to re-establish its equilibrium configura- 
tion. Note that a similar effect will occur in models of 
nonlinear capacitance. In that case the discharge time 
of the capacitor is an RC time constant, which clearly 
grows with external resistance R. 23 

The nonlincaritics of the storage time dependence on 
forward bias current and resistance have been noticed in 
the past. However, we are not aware of any group noting 
the following nonlinearity in the reverse recovery time. 
When the circuit in the inset of Fig. 2(b) is used with a 
repetitive square wave pulse sequence, one finds a strong 
dependence of the reverse recovery time on the pulse fre- 
quency. We consider the case where the pulses remain 
sufficiently far apart that the diode appears to fully re- 
cover from the previous pulse before the next arrives. As 
the rate of pulsing increases, the recovery time dramati- 
cally decreases. The effect is illustrated in Figs. 5 and 
6. Fig. 6(b) shows recovery traces from an NTE610 
diode at two pulse frequencies. With a pulse repetition 
rate of 200 kHz, the storage time is about 40 ns, while at 
a repetition rate of 3 MHz it has decreased to less than 
20 ns. Figs. 5 and 6 show that in general the storage 
time is constant at low frequencies, and then begins to 
fall roughly like f -1 / 2 with increasing frequency. This be- 
havior is seen in both experiment and PSpice simulation 
of the 1N4002, NTE610, and Linsay 1 model diodes. 
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FIG. 5. Frequency dependence of storage time (for diodes) 
or discharging effect (for PLC cases). Here we use a square 
wave of amplitude 1.5 V (3 V peak-to-peak), no DC bias, 50 
% duty cycle, (a) Diode 1N4007 and R=25 SI. Solid line 
is the experimental result, dashed line is the 1N4002 PSpice 
model, and dashed-dotted line is a 1/f 1 ^ 2 dependence, (b) 
PLC SPICE simulations of discharging storage time using pa- 
rameters from Ref. 4. The nonlinear resistor (NLR) uses the 
I-V curve of Ref. 3, while the PLC+R uses a R = 10 SI par- 
allel resistor. The Linsay model is a PSpice simulation using 
the diode model of Ref. 1 with R = 25 SI. Also shown as 
dashed lines are 1/f and 1/f 1 / 2 dependencies. 

From examination of charge dynamics in the p-n junc- 
tion it is known that the diffusion capacitance Co is a 
function of frequency w for sinusoidal driving signals; 9 



C D 




(4) 



where g is the differential conductance of the diode g 
= dl/dV. The f -1 / 2 frequency dependence of the capaci- 
tance comes from the diffusive nature of charge transport 
in the junction region. 9 The diffusion capacitance dimin- 
ishes with increasing frequency beyond f ~ l/27rr m , re- 
sulting in less charge stored in the p-n junction for a given 
forward bias current. This in turn implies that less mi- 
nority carrier charge must be moved across the junction 
at high frequencies, resulting in a shorter storage and re- 
verse recovery time. This is consistent with the results 
presented in Figures 5 and 6 for two rather different 
diodes (1N4007 and NTE610). 

One possible consequence of the frequency-dependent 
storage time is the persistence of period doubling behav- 
ior to unexpectedly high frequencies in the driven RLD 
circuit. As the driving frequency increases, the diffu- 
sion capacitance and storage time decrease, effectively 
increasing both the resonant frequency fo and 1/trr. 
This nonlinearity essentially helps to maintain the condi- 
tion f <~ fo ~ 1/ trr (where period doubling and chaos is 
most easily observed) over an extended range of driving 
frequencies. Thus period doubling and chaotic behavior 
should be more robust in the driven RLD circuit than 
one might naively expect. 

Fig. 6(a) illustrates the duty-cycle dependence of the 
storage time and its frequency dependence. The duty 
cycle is defined as the ratio of the time spent at the pos- 
itive amplitude voltage value to the total period of the 
periodic square pulse signal. As the duty cycle increases, 
there is greater overlap between adjacent pulses, causing 
the enhancement of the reverse recovery time, as origi- 
nally noted by Hunt and collaborators and summarized 
in Eqs. (l)-(3). Again we see very consistent behavior in 
both experiment and PSpice simulations. 
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FIG. 6. (a) Frequency and duty cycle dependence (20%, 
50%, and 80%) of the storage time for a NTE610 diode 
with R=25 Q, a square wave amplitude of 1.5 V (3 V 
peak-to-peak), and no DC bias. Solid lines are experimen- 
tal results and dashed lines are for a calculation using the 
FMMV2101 PSpice Zetex model, (b) Experimental time se- 
ries of diode current for the 50% duty cycle square wave case 
(from item (a)) for frequencies of 200 kHz and 3 MHz, illus- 
trating the decrease in reverse recovery time as the frequency 
increases. 

Finally, we have noted a strong DC voltage bias depen- 
dence of the storage time. The data and simulations on 
two different diodes is shown in Fig. 7. For forward bias 
the storage time increases dramatically, while for nega- 
tive bias it decreases. The range of storage time variation 
is more than a factor of 10 for a 2 volt swing in DC bias 
voltage. A forward-biased diode will have a broader re- 
gion of minority charge distribution, forcing the carriers 
to travel greater distances. This in turn increases the 
recovery time of the diode. Similarly, a reverse biased 
diode is able to sweep the minority carriers across the 
junction more quickly with the additional driving force 
provided by the reverse DC bias. This strong nonlinear- 
ity of the recovery time has not been discussed in the 
nonlinear dynamics literature, to our knowledge. 
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FIG. 7. DC voltage offset dependence of the storage 
time for two diodes. R=25 Q, amplitude of 1.5 V (3 V 
peak-to-peak), f=20 kHz, 50 % duty cycle, (a) Solid line is 
the experimental result on NTE610 diode and the dashed line 
is for a PSpice simulation using the FMMV2101 Zetex model, 
(b) Solid line is for experimental result on 1N4007 diode and 
dashed line is for PSpice simulation using FMMV2101 Zetex 
model. 

The DC bias nonlincarity can have a significant effect 
on the bifurcation diagram of the driven RLD circuit. 24 ' 25 
As a DC bias is applied, the storage time and differen- 
tial capacitance change, and the relationship between f , 
and the driving frequency f will change. As dis- 
cussed above, this can influence whether or not chaos is 
observed in the circuit. Here we present an example of 
the complex dynamics that can result. Figure 8 shows a 
case where the driving frequency f is larger than 1/trr 
and about two to three times the RLD circuit resonant 
frequency fo at zero DC bias. In this bifurcation diagram 
we sec only period- 1 behavior for driving signal power up 
to +15 dBm for both zero and positive DC bias. How- 
ever, for negative bias, we see period doubling set in at 
about +14 dBm for a -0.4 V bias and about +11 dBm for 
a -0.6 V bias. The prevailing NLC models cannot explain 
this observation because a negative bias moves the diode 
away from the nonlinear part of the C(V) curve. In this 
case we believe that the increase in 1 /trr and increase in 
fo due to Vdc create a situation where f ~ f ~ 1/trr, 
and the circuit displays period doubling at frequencies 
where ordinarily it only shows linear response. 
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FIG. 8. Bifurcation diagrams for RLD circuit (Fig. 1(a)) 
using an NTE610 diode with R=25 SI, f=29 MHz, and L=10 
/iH (note that f is greater than f = 12.3 MHz), (a) DC Offset 
= +0.2 V, (b) DC Offset = V, (c) DC Offset = -0.4 V, (d) 
DC Offset = -0.6 V. 

It should be noted that period doubling and chaos can 
be observed in models of the driven RLD circuit even in 
the apparent absence of reverse recovery. 3,10,26 However, 
it is not clear to what extent these models may have 
discharge-time effects that mimic the behavior of reverse 
recovery in diodes. 



simply reproduce the essential behavior.) 

The history dependence of the discharge time for the 
PLC model is illustrated in Fig. 3(b) (all simulations of 
the PLC model are done in SPICE). The discharge time 
is clearly enhanced after the second pulse of a two pulse 
sequence, compared to just a single pulse. This result 
is similar to that in Ref. 18 and Fig. 3(a) on real diodes. 
There is also strong frequency dependence to the dis- 
charge time of the PLC model, as illustrated in Fig. 5(b). 
The PLC capacitor model alone shows a discharge time 
that varies as 1 /f . When a linear resistor is added in par- 
allel with the PLC, the discharge time drops more slowly 
with frequency (PLC +R in Fig. 5(b)). Finally, when a 
nonlinear resistor (NLR) is added in parallel to the PLC, 
using the ideal diode IV curve of Linsay 1 , the discharge 
time shows behavior remarkably similar to the reverse 
recovery time of a real diode. 

In summary, the discharge behavior of the simple PLC 
model (when implemented in the circuits of interest) is 
remarkably similar to the reverse recovery behavior of 
real diodes. We find that SPICE simulations of other 
more sophisticated NLC models give very similar results. 
We conclude that all NLC models will show a reverse- 
recovery-like effect. It appears that the discharge of the 
NLC/PLC model is essentially the same as the diffusion 
of charge carriers out of the p-n junction. Hence this 
important ingredient of nonlinear physics is common to 
the two major models of nonlinear dynamics of the driven 
nonlinear resonator. However, our results demonstrate 
that all of these models lack key additional nonlinear 
properties of the reverse-recovery time. 



VI. CONCLUSIONS 



V. DISCUSSION 

To make further comparisons between the nonlinear 
capacitor (NLC) and reverse-recovery time models of 
RLD dynamics, we wish to further consider the simpli- 
fied PLC model of the NLC. The PLC model contains 
the essential physics of all NLC models, and is very easy 
to implement and understand. Figure 2(a) (dashed line) 
shows that the PLC model of the diode (Fig. 3(b), 
inset) clearly shows an effective delay in "turning off' 
after a large forward current has been applied to the de- 
vice. This delay is a combination of two RC time con- 
stants. The first, or "storage" part of the response is 
associated with the RC discharge of the PLC during its 
high-capacitance phase. As the voltage drop on the PLC 
decreases, it reaches the point (V = E ) where the capac- 
itance switches to its smaller value, resulting in a much 
more rapid discharge, resembling the transition phase of 
reverse recovery. The simulation in Fig. 2(a) (dashed 
line) for the PLC qualitatively resembles the data on the 
1N4007 diode. (Note that this simulation is not intended 
to fit the precise behavior of the 1N4007 diode, but to 



We maintain that reverse recovery is a key effect pro- 
ducing chaos in the sinusoidally driven RLD circuit. All 
nonlinear capacitor models of this circuit have a behavior 
that mimics reverse recovery. Both the reverse-recovery 
and nonlinear capacitor models have memory built into 
them through a charge storage mechanism. The details 
of these models are different, but the results for circuit 
dynamics are basically the same. We find that it is not 
useful to make the distinction that either the nonlinear 
capacitance or the reverse recovery on minority carriers is 
the cause of chaos, because they can have similar effects 
on the circuit dynamics. Finally, the reverse recovery 
time in real diodes is a strongly nonlinear function of 
forward bias current, driving frequency and DC bias. A 
detailed understanding of chaos in the RLD circuit must 
include the reverse recovery effect and all of its nonlin- 
earities. 
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TABLE I. Survey of results for period doubling and chaos on driven RLD circuits with four different diodes. The reverse 
recovery time trr was measured at 20 kHz (where it is frequency independent) with V dc = 0, a 1.5 V amplitude pulse (3 V 
peak-to-peak), 50 % duty cycle, and R = 25 Q. The zero-bias junction capacitance Cj was measured with an Agilent 4285A 
LCR meter. The resonant frequencies fo were changed by adjusting the L value in the RLD circuit, while R was fixed at 25 Q. 
The search for period doubling and chaos took place over the range f ~ 0.1 fo to f ~ 10 fo. We were unable to measure a finite 
reverse recovery time (i.e. trr < 10 ns) in the 1N4148 fast-recovery diode due to the bandwidth limitations of our oscilloscope. 
Thus we were unable to establish the condition f ~ f ~ 1/trr for that diode. 





Diods 


trr l,nsj 


(p? ) 


lU-bLUlb Willi 


IXcbLllLb Willi 


llCD III Lb Wl 111 








fo ~ 1/trr 


fo ~ W/trr 


fo ~ WO/trr 


1N5400 


7000 


81 


Period doubling and 


Period doubling and 


No period doubling 








chaos for 


chaos for 


or chaos 








f/fo ~ 0.11 - 1.64 


f/fo ~ 0.16 - 1.76 




1N4007 


700 


19 


Period doubling and 


Period doubling and 


No period doubling 








chaos for 


chaos for 


or chaos 








f/fo ~ 0.13 - 2.0 


f/fo ~ 0.23 - 1.3 




1N5475B 


160 


82 


Period doubling and 


No period doubling 


No period doubling 








chaos for 


or chaos 


or chaos 








f/fo ~ 0.66 - 2.2 






NTE610 


45 


16 


Period doubling and 


Period doubling 


No period doubling 








chaos for 


only for 


or chaos 








f/fo ~ 0.14 - 3.84 


f/fo ~ 1.17 - 3.25 
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